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SECTION  I 

INTRODUCTION  AND  SUMMARY 


This  report  discusses  the  results  using  various  £-k  spectral 
techniques  in  locating  live  fairly  weak  events  with  magnitudes  ranging  from 
4.4  to  4.8.  Three  different  subarray  configurations  (Figure  I-l)  were  used 

to  process  these  events: 

•  Processed  subarray  outputs  from  AO  through 
the  D-ring 

•  The  extended  E3  subarray 

•  The  standard  Fl  subarray 


To  complement  theoretical  f-k  spectral  work  using  location 
statistics,  '  this  experiment  applied  various  f-k  spectral  techniques  to  some 
of  the  seismic  events.  There  are  several  pertinent  questions  to  be  answere 


d. 


•  Does  any  £-k  technique  give  a  more  accurate 
epicentral  location  than  the  others? 

•  Can  some  £-k  techniques  give  accurate  location 
while  requiring  significantly  fewer  calculations 
than  other  techniques? 

•  How  do  the  epicentral  locations  obtained  from 
arrays  of  different  apertures  compare? 

It  was  found  that  the  best  strategy  for  approximately  deter¬ 
mining  the  epicenter  of  weak  events  at  LASA  is  to  use  bandpassed  subarray 
outputs  to  form  multiple  beams  in  the  time  domain.  The  extended  E3  sub¬ 
array  gave  correct  approximate  epicenters  for  the  events  studied,  while  a 
standard  subarray  does  not  give  satisfactory  approximate  epicenters. 


;  T  ?  T  ?-r- 

LASA 


Figure  I-l.  Array  Configurations 
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A.  SHORT-PERIOD  EVENTS 

Five  short-period  events  of  magnitudes  ranging  from  4.4  to 
d.8  were  selected  to  be  processed.  Table  I- 1  describes  the  selected  events. 
Even,  locations  were  determined  from  PDF  information  supplied  by  the  Coa, 
and  Geodetic  Survey  and  the  LASA  event  log.  Two  of  the  f.ve  events  com 
not  be  identified;  their  epicenter  locations  were  estimated  from  the  LAS 

subarray  outputs. 

B.  ARRAY  DESCRIPTION 

The  larges,  subarray  con.-isted  of  outputs  from  the  LASA  AO, 

D-ring  subarray.  Outer  LASA  E-  and  F-rings  were  no,  included  in  frequency- 

wavenumber  analysis  because  the.r  traveitime  anomalies  and  long  moveou, 

times  badly  disrupt  the  frequency-wavenumber  analysis.  Tbe  intermed.ate 

subarray  was  composed  of  the  extended  E3  subarray,  and  the  smallest  sub- 

of  the  FI  subarray.  In  the  latter  cases,  all  25  seis- 
array  was  composed  ot  tne  r  i  suunn  y 

mometer  outputs  were  used. 

C.  EVENT  DETECTION 

Even,  location  presuppose,  detection.  E.ven,  detection  was  ac¬ 
complished  by  digitally  simulating  the  LASA  on-line  telesei.m  detector  using 
the  subarray  sums  which  had  been  prefiitered  to  eliminate  the  microsei.m 
peak  and  pass  the  frequency  band  generally  containing  most  of  the  P  wave 

energy. 

An  independent  event  detector  was  connected  with  each  subarray 
output  (sensors  were  connected  with  Ei.  E2,  E4,  FI,  F2.  F3.  and  F4).  Sin¬ 
gle-channel  detection  consi,  ted  of  rectifying  the  channel  output  and  passing  the 
rectified  waveform  through  a  low-pass  recursive  filter  of  the  form 


V  .  ,  =  a  (I  -  a)  A  +  (1  -  a)  X 

n  +  I  " 
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SHORT-PERIOD  PDE  INFORMATION 


where 


-1/300 

a  =  e 

A  ,  =  filter  output  at  the  nrl*^^  point 
n+  1 

X  =  point  of  single-channel  time  series 

n 


The  most  recent  rectified  time  sample,  X^,  then  was  compared 
with  the  recursive  filter  output  delayed  by  60  sec  (i.  e.  ,  600  samples).  If  this 
ratio  is  greater  than  5.  82,  the  decision  is  made  than  an  event  has  occurred; 


i.  e.  , 

X  /{A  -  600)  >5,82 

n  /  n 


If  four  out  of  the  seven  (El,  E2,  E4,  FI,  F2,  F3,  and  F4)  single-channel  de¬ 
tectors  indicate  that  an  event  is  present  witlan  a  20-sec  interval,  the  overall 
decision  is  that  a  teleseism  has  occurred. 


The  results  of  applying  this  detection  technique  are  shown  in 
Figures  1-2  through  1-6.  All  21  subarray  outputs  are  shown  after  filtering 
but  only  the  E-  and  F-ring  (less  E3)  subarray  outputs  were  used  in  the  process. 

Individual  trace  detectin'^  intervals  are  shown  as  a  first  point 
of  detection  in  the  trace  playbacks.  Detection  points  are  indicated  by  an  ar¬ 
row;  however,  traces  from  E  and  F  subarrays  which  did  not  trigger  a  detec¬ 
tor  are  not  marked.  All  five  teleseisrns  were  successfi’lly  detected  using  the 
four-out-of-seven  decision  rule. 
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Figure  1-2.  Filtered  Subarray  Output*,  Event  1 
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Figure  1-3,  F-Ring  Subarray  Outputs,  Lvent  2 
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Figure  1-5.  F-Ring  Subarray  Outputs,  Event  4 
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Figure  1-6.  F-Rlng  Subarray  Outputs,  Event  5 
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SECTION  II 


MATHEMATICAL  DESCRIPTION  OF  f-k  SPECTRA 


The  conventional  spectrum  is  defined  as 
P^(V)=  V  FF  V 


where 

V  is  the  vector  of  phase  shifts  appropriate  for  the 
frequency  and  velocity  of  a  plane  wave  propa¬ 
gating  across  the  array 

* 

denotes  conjugate  transpose 

F  is  the  vector  of  Fourier  transforms  at  the  desired 
frequency  at  each  sensor  which  are  taken  in  the 
same  order  as  ^ 


These  transforms  are  taken  over  a  data  gate  including  the 
P-wave  signal  and  are  normalized  (only  phase  information  is  retained)  to 
remove  amplitude  anomalies.  This  may  be  interpreted  as  a  beamsteer  cal¬ 
culated  in  the  frequency  domain. 


P2{V)  is  defined  as 


V  f  f  V 


where  f  is  a  vector  of  filter  weights  determined  by  the  equation 


j^kl  +  FF*]  |f|  = 


0 

0 

k 

0 

0 

0 


(2-1) 


II- 1 
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where 


k  is  a  scaler  (usually  0.  01) 
I  Is  the  Identity  matrix 


the  position  of  k  in  the  right-hand  side  of  equation 
2-1  indicates  which  channel  is  used  as  reference 


This  process  can  be  thought  of  as  designing  a  spatial  whitening  filter  and  then 
considering  the  reciprocal  of  the  power  response  in  wavenumber  space  of  the 
multichannel  filter  system. 

P^(V)  is  defined  as 


k 


i=l 


1 


i  1 

where  the  f.  's  are  the  filter  weights  designed  with  the  i^^sensor  used  as  ref 
erence  in  Equation  2-1.  This  spectral  estimate  is  obtained  using  several 
sensors  as  reference.  If  all  sensors  are  used  as  reference,  the  result  is 

^  V  (kl  +  FF*)  V 


Pg(V)  is  defined  as 


1 


V  (kl  +  FF  )  V 

which  can  be  obtained  in  the  following  manner.  Consider  D  =  N  -  MW  ,  where 
N  is  a  measured  covariance  matrix  (kl  +  FI  )  and  M  is  a  scalar.  The  prob¬ 
lem  is  to  find  the  largest  M  such  that  D  will  remain  a  positive  definite  matrix. 
D  remains  ii  positive  definite  matrix  for 


M  = 


-1  -* 
V  N  V 


U-2 
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This  is  one  measure  of  how  large  a  portion  of  the  matrix  N  consists  of  mode 
V^,  The  expression 


-♦*  -1 
V  N  V 


where  N'  ^  =  (kl  +  FF  *)  \  also  can  be  thought  of  as  the  output  of  a  maximum- 


likelihood  filter  designed  to  pass  the  plane-wave  signal 
having  a  covariance  matrix  of  (kl  -f-  FF  ). 


in  a  noise  field 


It  should  be  noted  that  for  the  preceding  definitions  and  for  us¬ 
ing  N  =  (kl  -f  FF*),  Pj(V),  P4(V),  and  Pg(V)  will  necessarily  have  peaks  at 
the  same  V  (Appendix  A).  Since  this  study  was  principally  interested  in  the 

H  H  -4 

peaks  for  location  purposes,  only  estimates  of  Pj^(V),  P2(V),  and  P^(V)  were 
computed. 


n-3/4 
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SECTION  m 
PROCESSING  RESULTS 


Five  short-period  events  were  processed  using  a  15-sec  gate 
for  the  LASA  AO,  D-ring  subarray,  a  10-sec  gate  for  the  extended  E3  sub¬ 
array,  and  a  5-sec  gate  for  the  FI  subarray.  Gate  position  locations  were 
based  on  time-trace  playbacks  and  were  chosen  to  include  the  primary  P-wave 
arrivals. 

Plots  for  all  events  were  made  at  fixed  frequencies  of  0.  850, 

1.  000,  and  1,  115  Hz,  An  additional  plot  was  made  at  peak  transform  frequen¬ 
cy  computed  from  the  center  sei-smometei*  time  gate. 

The  subarray  outputs  used  to  compute  the  AO,  D-ring  f-k 
spectra  were  a  summation  of  the  center  element  and  rings  2,  3,  4,  and  5 
(13  elements).  This  configuration  was  chosen  because  of  good  isotropic  sig¬ 
nal  response  (V  >10  km/sec)  and  fairly  good  noise  rejection  capabilities  in 
the  0.  6  Hz  <  f  <  1.  5  Hz  range. 

For  the  E3  and  FI  subarray  data,  the  inputs  were  the  25  indi¬ 
vidual  sensors.  In  every  case  the  data  was  low-cut  filtered  before  taking  the 
transforms. 

Results  of  the  frequency-wavenumber  processing  analysis  are 
summarized  in  Tables  B-1  through  B-11  in  Appendix  B.  These  tables  give 
the  location  of  the  f-k  spectrum's  peak  for  each  type  of  estimate  at  each  fre¬ 
quency.  The  tables  also  contain  comments  about  the  behavior  of  the  spectral 
estimate  in  the  f-k  region  of  the  incident  teleseism.  Peak  locations  and  epi¬ 
centers  are  specified  by  dividing  the  frequency -wave number  plots  into  a 
3000-point  grid  numbered  from  the  lower  left-hand  corner.  There  are  49  grid 


IU-1 


8ci«noe  ••rvicas  division 


values  in  the  k  direction  and  59  grid  values  in  the  k  direction,  with  the  plot 

X  y 

origin  at  the  center  of  the  grid.  Consequently,  an  event  arriving  with  a  ver¬ 
tical  incidence  would  be  located  at  k  =  25,  k  =  30. 

X  y 

All  f-k  spectra  calculated  in  this  experiment  are  not  reproduced 
in  this  report.  Some  spectra  from  the  event  recorded  on  21  June  1967  from 
eastern  Russia  or  northeastern  China  are  shown  as  examples.  Figure  lU-l 
shows  the  three  types  of  f-k  spectra  calculated  at  0.  85  He  (the  peak  power  of 
the  center  array  transform).  All  f-k  spectra  give  a  reasonable  accurate  lo¬ 
cation  at  the  "peak"  frequency  for  this  event;  however,  this  is  not  true  for  all 
events. 

•4 

Figure  III-2  shows  the  three  types  of  f-k  spectra  calculated  at 
1.  00  He.  At  this  frequency  none  of  the  f-k  spectra  peak  at  the  proper  k,  but 
conventional  and  high-resolution  spectra,  using  all  the  outer  ring  (D-ring) 
subarrays  as  reference,  did  give  some  indication  of  the  event. 
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Figure  III-l.  Frequency-Wavenumber  Plots  of  21  June  1967  Event  at  0.  85  Hz 
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SECTION  iV 

DISCUSSION  OF  ON-LINE  EVENT  LOCATION 

This  experiment  was  intended  to  indicate  whether  or  not  f-k 
spectral  techniques  offer  any  advanUges  for  approximate  on-line  location. 

Such  an  on-line  location  scheme  must  have  a  predetermined  mode  of  opera¬ 
tion  which  is  either  entirely  fixed  or  dependent  only  upon  relatively  simple 
calculations  based  on  incoming  data. 

Assuming  detection,  three  methods  of  operation  seem  reasonable. 

•  Calculating  f-k  spectra  at  a  few  fixed  frequen¬ 
cies  which  s^jan  the  band  of  most  energy. 

•  Calculating  peak  energy  frequency  in  the  signal 
gate  and  then  calculating  an  f-k  spectrum  at  this 
frequency 

•  Calculating  f-k  spectra  over  a  range  of  frequen¬ 
cies  and  executing  a  velocity-preserving  stack 
of  the  resulting  spectra 

In  the  third  method,  this  process  can  be  done  more  efficiently 
in  the  time  domain  for  the  conventional  spectrum  as  the  number  of  frequen¬ 
cies  in  the  stack  increases.  In  the  time  domain,  this  method  becomes  a 

*4 

bandpass -filter  and  form-many-beams  operation.  An  f-k  approach  in  this 
case  might  be  advantageous  only  if  one  of  the  high-resolution  techniques  proved 
significantly  superior. 

The  best  approximate  location  at  LASA  clearly  can  be  obtained 
using  the  large  array  composed  of  subarray  outputs.  These  outputs  have  two 
advantages  over  individual  subarrays;  greater  aperature  and  better  S/N  ratios 
(on  the  average).  Data  using  the  large  array  from  AO  through  the  D-ring  sug¬ 
gest  that  f-k  epicenter  location  could  be  done  most  reliably  by  forming  velocity- 
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preierving  itacki.  Th«  data,  further  indicate  that  the  conventional  P^(V)  spec¬ 
tra  are  genernllv  more  reliable  than  high-reeolution  spectra  which  are  the 
power  response  of  spatial  whitening  filters  [P2(V)or  Pj(V)].  The  implica¬ 
tion  of  these  data  are  that  the  best  strategy  for  roughly  locating  events  is  to 
use  bandpassed  data  a.nd  form  beams  in  the  time  domain  to  give  a  power  den¬ 
sity  picture  in  the  wavenumber  or  velocity  plane. 

Data  from  the  extended  E3  subarray  were  limited  to  three  events. 
In  these  cases  the  events  could  have  been  located  approximately  by  using  either 
the  conventional  spectrum  at  the  transform's  peak  frequency  or  a  velocity  pre- 
serving  stack.  These  data  also  indicate  that  the  high-resolution  f-k  spectra 
[P^CV)  +  o^®r  poorer  event  location  than  do  the  conventional 

spectra  (Pj(V)). 

Data  from  the  FI  subarray  were  limited  to  three  events.  In  all 
three  cases  no  reasonable  approximate  epicenter  could  be  obtained  from  the 
FI  subarray  data. 
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APPENDIX  A 


LOCATION  OF  MAXIlVlA  FOR  P^(V),  P^CV),  and  Pg(V) 

The  following  values  are  given  for  Pj(V),  P^(V),  and  Pg(V): 


P^{V)  =  V  FF  V 
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P  .(V)  =  - 

4  -»4i 
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-2"* 

V  (kl  +  FF  )  V 
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(kl  +  FF*)'^  and  (kl  +  FF*)'^  naay  be  stated  as 
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+  J.  fl 
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(k 
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Substituting  these  values  in  Equations  A-2  and  A-3  yields 


P^(V) 


V  V  - 
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(k+  F  F) 


2kV  FF  V  + 
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F(F  F)F  V 
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{A-2) 


Pc(V)  = 


V  V  - 


k  + 


[V  FF  V] 
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Since 

V  (kl  +  FF  )  V  >  0 

-2  *♦ 

V  (kl  +  FF  )  V  >  0 

and 

V  V  K  a  constant  (the  number  of  channels) 

■4  e-f 

F  F  s  a  constant  (independent  of  V) 

it  follows  that  Pj(V),  and  Pg(^)  will  all  have  a  maximum  at  the  maxi¬ 

mum  of  (V*F?*^). 
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